Genetic variation among 82 pharmacogenes: The PGRNseq data from the eMERGE network.

Genetic variation can affect drug response in multiple ways, although it remains unclear how rare genetic variants affect drug response. The electronic Medical Records and Genomics (eMERGE) Network, collaborating with the Pharmacogenomics Research Network, began eMERGE-PGx, a targeted sequencing study to assess genetic variation in 82 pharmacogenes critical for implementation of "precision medicine." The February 2015 eMERGE-PGx data release includes sequence-derived data from ∼5,000 clinical subjects. We present the variant frequency spectrum categorized by variant type, ancestry, and predicted function. We found 95.12% of genes have variants with a scaled Combined Annotation-Dependent Depletion score above 20, and 96.19% of all samples had one or more Clinical Pharmacogenetics Implementation Consortium Level A actionable variants. These data highlight the distribution and scope of genetic variation in relevant pharmacogenes, identifying challenges associated with implementing clinical sequencing for drug treatment at a broader level, underscoring the importance for multifaceted research in the execution of precision medicine.

A genome-wide association study identifies variants in KCNIP4 associated with ACE inhibitor-induced cough.

The most common side effect of angiotensin-converting enzyme inhibitor (ACEi) drugs is cough. We conducted a genome-wide association study (GWAS) of ACEi-induced cough among 7080 subjects of diverse ancestries in the Electronic Medical Records and Genomics (eMERGE) network. Cases were subjects diagnosed with ACEi-induced cough. Controls were subjects with at least 6 months of ACEi use and no cough. A GWAS (1595 cases and 5485 controls) identified associations on chromosome 4 in an intron of KCNIP4. The strongest association was at rs145489027 (minor allele frequency=0.33, odds ratio (OR)=1.3 (95% confidence interval (CI): 1.2-1.4), P=1.0 × 10(-8)). Replication for six single-nucleotide polymorphisms (SNPs) in KCNIP4 was tested in a second eMERGE population (n=926) and in the Genetics of Diabetes Audit and Research in Tayside, Scotland (GoDARTS) cohort (n=4309). Replication was observed at rs7675300 (OR=1.32 (1.01-1.70), P=0.04) in eMERGE and at rs16870989 and rs1495509 (OR=1.15 (1.01-1.30), P=0.03 for both) in GoDARTS. The combined association at rs1495509 was significant (OR=1.23 (1.15-1.32), P=1.9 × 10(-9)). These results indicate that SNPs in KCNIP4 may modulate ACEi-induced cough risk.

Genetic variation in the HLA region is associated with susceptibility to herpes zoster.

Herpes zoster, commonly referred to as shingles, is caused by the varicella zoster virus (VZV). VZV initially manifests as chicken pox, most commonly in childhood, can remain asymptomatically latent in nerve tissues for many years and often re-emerges as shingles. Although reactivation may be related to immune suppression, aging and female sex, most inter-individual variability in re-emergence risk has not been explained to date. We performed a genome-wide association analyses in 22,981 participants (2280 shingles cases) from the electronic Medical Records and Genomics Network. Using Cox survival and logistic regression, we identified a genomic region in the combined and European ancestry groups that has an age of onset effect reaching genome-wide significance (P>1.0 × 10(-8)). This region tags the non-coding gene HCP5 (HLA Complex P5) in the major histocompatibility complex. This gene is an endogenous retrovirus and likely influences viral activity through regulatory functions. Variants in this genetic region are known to be associated with delay in development of AIDS in people infected by HIV. Our study provides further suggestion that this region may have a critical role in viral suppression and could potentially harbor a clinically actionable variant for the shingles vaccine.

Design and anticipated outcomes of the eMERGE-PGx project: a multicenter pilot for preemptive pharmacogenomics in electronic health record systems.

We describe here the design and initial implementation of the eMERGE-PGx project. eMERGE-PGx, a partnership of the Electronic Medical Records and Genomics Network and the Pharmacogenomics Research Network, has three objectives: (i) to deploy PGRNseq, a next-generation sequencing platform assessing sequence variation in 84 proposed pharmacogenes, in nearly 9,000 patients likely to be prescribed drugs of interest in a 1- to 3-year time frame across several clinical sites; (ii) to integrate well-established clinically validated pharmacogenetic genotypes into the electronic health record with associated clinical decision support and to assess process and clinical outcomes of implementation; and (iii) to develop a repository of pharmacogenetic variants of unknown significance linked to a repository of electronic health record-based clinical phenotype data for ongoing pharmacogenomics discovery. We describe site-specific project implementation and anticipated products, including genetic variant and phenotype data repositories, novel variant association studies, clinical decision support modules, clinical and process outcomes, approaches to managing incidental findings, and patient and clinician education methods.

Branched-chain amino acid levels are associated with improvement in insulin resistance with weight loss.

AIMS/HYPOTHESIS: Insulin resistance (IR) improves with weight loss, but this response is heterogeneous. We hypothesised that metabolomic profiling would identify biomarkers predicting changes in IR with weight loss. METHODS: Targeted mass spectrometry-based profiling of 60 metabolites, plus biochemical assays of NEFA, ß-hydroxybutyrate, ketones, insulin and glucose were performed in baseline and 6 month plasma samples from 500 participants who had lost ≥4 kg during Phase I of the Weight Loss Maintenance (WLM) trial. Homeostatic model assessment of insulin resistance (HOMA-IR) and change in HOMA-IR with weight loss (∆HOMA-IR) were calculated. Principal components analysis (PCA) and mixed models adjusted for race, sex, baseline weight, and amount of weight loss were used; findings were validated in an independent cohort of patients (n = 22). RESULTS: Mean weight loss was 8.67 ± 4.28 kg; mean ∆HOMA-IR was -0.80 ± 1.73, range -28.9 to 4.82). Baseline PCA-derived factor 3 (branched chain amino acids [BCAAs] and associated catabolites) correlated with baseline HOMA-IR (r = 0.50, p < 0.0001) and independently associated with ∆HOMA-IR (p < 0.0001). ∆HOMA-IR increased in a linear fashion with increasing baseline factor 3 quartiles. Amount of weight loss was only modestly correlated with ∆HOMA-IR (r = 0.24). These findings were validated in the independent cohort, with a factor composed of BCAAs and related metabolites predicting ∆HOMA-IR (p = 0.007). CONCLUSIONS/INTERPRETATION: A cluster of metabolites comprising BCAAs and related analytes predicts improvement in HOMA-IR independent of the amount of weight lost. These results may help identify individuals most likely to benefit from moderate weight loss and elucidate novel mechanisms of IR in obesity.

Genome-wide association scan for childhood caries implicates novel genes.

Dental caries is the most common chronic disease in children and a major public health concern due to its increasing incidence, serious health and social co-morbidities, and socio-demographic disparities in disease burden. We performed the first genome-wide association scan for dental caries to identify associated genetic loci and nominate candidate genes affecting tooth decay in 1305 US children ages 3-12 yrs. Affection status was defined as 1 or more primary teeth with evidence of decay based on intra-oral examination. No associations met strict criteria for genome-wide significance (p < 10E-7); however, several loci (ACTN2, MTR, and EDARADD, MPPED2, and LPO) with plausible biological roles in dental caries exhibited suggestive evidence for association. Analyses stratified by home fluoride level yielded additional suggestive loci, including TFIP11 in the low-fluoride group, and EPHA7 and ZMPSTE24 in the sufficient-fluoride group. Suggestive loci were tested but not significantly replicated in an independent sample (N = 1695, ages 2-7 yrs) after adjustment for multiple comparisons. This study reinforces the complexity of dental caries, suggesting that numerous loci, mostly having small effects, are involved in cariogenesis. Verification/replication of suggestive loci may highlight biological mechanisms and/or pathways leading to a fuller understanding of the genetic risks for dental caries.

Discovery, purification, and characterization of a temperate transducing bacteriophage for Bordetella avium.

We discovered and characterized a temperate transducing bacteriophage (Ba1) for the avian respiratory pathogen Bordetella avium. Ba1 was initially identified along with one other phage (Ba2) following screening of four strains of B. avium for lysogeny. Of the two phage, only Ba1 showed the ability to transduce via an allelic replacement mechanism and was studied further. With regard to host range, Ba1 grew on six of nine clinical isolates of B. avium but failed to grow on any tested strains of Bordetella bronchiseptica, Bordetella hinzii, Bordetella pertussis, or Bordetella parapertussis. Ba1 was purified by CsCl gradient centrifugation and was found to have an icosahedral head that contained a linear genome of approximately 46.5 kb (contour length) of double-stranded DNA and a contractile, sheathed tail. Ba1 readily lysogenized our laboratory B. avium strain (197N), and the prophage state was stable for at least 25 generations in the absence of external infection. DNA hybridization studies indicated the prophage was integrated at a preferred site on both the host and phage replicons. Ba1 transduced five distinctly different insertion mutations, suggesting that transduction was generalized. Transduction frequencies ranged from approximately 2 x 10(-7) to 1 x 10(-8) transductants/PFU depending upon the marker being transduced. UV irradiation of transducing lysates markedly improved transduction frequency and reduced the number of transductants that were lysogenized during the transduction process. Ba1 may prove to be a useful genetic tool for studying B. avium virulence factors.